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Ri jke 管 热 声 不 稳定 非 线 性 动力 学 研究 
-分 岔 分 析 


zh, HX, AMA 
(西北 工业 大 学 ， 燃 烧 、 流 动 和 热 结 构 国 家 级 重点 实验 室 ， 西 安 710072) 

摘 要 : 建立 了 水 平 Rijke 管 热 声 模型 ， 利 用 Galerkin 方法 对 控制 方程 组 进行 展开 ， 实 现 数值 求解 。 确 定 了 Galerkin 
模 态 收敛 阶 数 为 10， 利 用 非 线 性 动力 学 理论 对 系统 进行 分 命 分 析 。Rijke 管 热 声 系 统 分 命 行为 属于 亚 临 界 Hopf 分 贫 。 
系统 稳定 性 区 域 分 为 全 局 稳定 、 全 局 不 稳定 及 双 稳 态 区 域 。 获 得 了 无 量 纲 加 热 功率 K, MARE yo IMERI cl 
和 时 间 延 迟 zt 等 参数 的 分 命 图 谱 ， 发 现 加 热 器 位 置 xy 的 分 命 图 谱 存 在 两 个 Hopf 分 岔 点 。 在 线性 不 稳定 区 域内 ， 振 
WEL BELEN TA) REIS t 的 增 大 呈现 先 增 大 后 减 小 的 趋势 。 
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Rijke Tube Thermoacoustic Instability Nonlinear 
Dynamics Study—Bifurcation Analysis 


FENG Jian-Chang AO Wen LIU Pei-Jin 
(Science and Technology on Combustion, Internal Flow and Thermal-structure Laboratory, 


Northwestern Polytechnical University, Xi’an 710072, China) 


Abstract: The thermoacoustic model of a horizontal Rijke tube is established. The governing equations 
are expanded and solved by using the Galerkin method. The convergence order of Galerkin modes is 
determined as 10. The bifurcation analysis of the system is carried out by using nonlinear dynamics theory. 
The bifurcation behaviors of the horizontal Rijke tube thermoacoustic system is found to be subcritical Hopf 
bifurcation. The system stability regions are divided into globally stable region, globally unstable region and 
bistable region. The bifurcation diagrams of the non-dimensional heater power (K), the position of the heater 
(xy), the damping coefficient (cı) and the time delay (t) are obtained, and the bifurcation diagram of the 
position of the heater (xp) has two Hopf bifurcation points. In linear unstable region, the amplitude of the 
oscillation firstly increases and then decreases with the increase of time lag (7). 
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0 引言 


燃烧 不 稳定 (combustion instability) 是 指 燃烧 
器 中 压力 和 速度 的 大 幅度 自 激 振荡 ， 这 一 振荡 带 来 
的 循环 负荷 会 造成 结构 过 度 振动 和 燃烧 室 壁 面 的 热 
量 过 度 传递 等 问题 ， 对 结构 组 件 造成 致命 的 破坏 。 
长 久 以 来 ， 燃 烧 不 稳定 一 直 都 是 火箭 发 动机 ， 燃 气 
轮机 ， 航 空 发 动机 和 电站 锅炉 等 许多 高 性 能 燃烧 装 
置 面临 的 重大 挑战 。 热 声 不 稳定 〈thermoacoustic 
instability) 是 燃烧 不 稳定 的 最 为 主要 的 一 种 形式 ， 
一 般 认 为 ， 热 声 不 稳定 是 由 于 不 稳定 放 热 和 声场 振 
荡 的 耦合 引起 的 ， 将 引发 接近 固有 声场 频率 的 大 振 
幅 压 力 振荡 ， 其 本 质 是 热源 的 放 热 与 声场 环境 存在 
ERHALE. 

燃烧 不 稳定 现象 背后 所 蕴含 的 机 理 十 分 复杂 ， 
包含 人 燃烧， 流动 和 声场 的 相互 耦合 作用 。 火 焰 面 的 
变化 ， 当 量 比 波动 ， 流 体 动力 学 不 稳定 导致 的 涡 脱 
落 以 及 燃料 的 雾 化 及 蒸发 过 程 中 的 振荡 都 可 能 驱动 
燃烧 不 稳定 。 而 燃烧 室 几 何 结构 , 燃料 /空气 混合 度 ， 
燃料 的 类 型 和 构成 , 以 及 工作 环境 (注入 空气 温度 ， 
燃烧 室 压 强 ， 当 量 比 ， 涡 流 强 度 ， 贫 油 预 混 等 ) 都 
将 对 燃烧 不 稳定 产生 影响 口 。 

流动 和 火焰 的 非 线性 动力 学 特性 是 燃烧 不 稳 
定 中 最 为 重要 的 现象 之 一 。 燃 烧 室 中 的 流动 和 火焰 
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用 导致 的 不 稳定 机 理 及 控制 的 基础 ， 因 此 有 许多 而 
究 是 针对 Rijke 管 来 开展 的 5658。 不 稳定 放 热 作为 热 
声 不 稳定 的 激励 ， 是 不 稳定 发 生 的 一 个 关键 环节 。 
热 声 系统 中 的 能 量 输入 可 以 有 多 种 形式 ， 例 如 燃烧 
室 中 放 热 的 化 学 反应 , 电热 丝 , 电热 薄膜 (hot gauze) 
“611, Jahnke 和 Culick0l0 最 早 将 现代 动力 系统 理论 
引入 燃烧 不 稳定 研究 之 中 ， 以 分 析 非 线性 燃烧 不 稳 
定 。 他 们 提出 了 一 种 延 拓 算法 ， 可 对 各 稳 态 和 极限 
环 的 解 进行 体系 化 高 效 计算 ， 由 此 可 得 到 稳 态 和 极 
限 环 下 的 稳定 性 ， 通 过 分 岔 分 析 可 以 判定 不 稳定 发 
生 点 。 不 稳定 发 生 点 处 的 声波 的 性 质 取 雇 于 分 岔 的 
类 型 和 分 命 点 处 极限 环 的 性 质 。,Ananthkrishnan 等 上 
针对 燃烧 室 中 的 不 稳定 运动 建立 了 降 阶 模型 ， 并 求 
解 分 析 了 触发 和 极限 环 性 质 。 分 析 表 明 ， 对 于 轴 向 
振荡 模 态 ， 一 阶 不 稳定 需要 至 少 四 阶 Galerkin 级 数 
展开 ， 而 二 阶 不 稳定 则 至 少 需 要 八 阶级 数 展开 ， 才 
能 确保 计算 的 准确 性 。Balasubramanian02 研 究 了 非 
正则 性 在 简单 热 声 系统 中 的 作用 和 影响 ， 通 过 推导 
水 平 Rijke 管 的 声场 控制 方程 并 求解 ， 重 点 讨论 非 
正则 性 和 非 线 性 效应 带 来 的 结果 。Subramanian[!3 
分 析 了 热 声 系 统 中 亚 临 界 分 岔 的 性 质 ， 以 及 随后 导 
致 的 双 稳 定性 问题 。Waugh 等 0445 研 究 了 热 声 系统 
中 的 触发 和 流体 力学 中 的 向 消 流 的 旁 路 转换 (逾越 
型 ”转换 ) 的 比拟 ， 二 者 机 理 均 是 小 扰动 导致 系统 


特性 随 着 控制 参数 的 变化 会 发 生 剧烈 改变 ， 尤 其 是 
当 参 数 变化 经 过 分 个 点 时 。 燃 烧 过 程 本 身 可 能 会 ， 
抑或 不 会 出 现 分 岔 现象 ， 但 一 旦 发 生 分 岔 ， 非 线性 
行为 在 许多 燃烧 装置 中 将 占据 主导 地 位 BH。 系统 参 
数 的 微小 扰动 有 可 能 会 引起 分 贫 ， 导 致 系统 由 稳 态 
〈 小 振荡 或 无 振荡 ) 转换 为 不 稳定 状态 〈 大 振幅 的 
极限 环 )。 相反 地 , 在 相同 的 临界 参数 值 时 ， 系 统 由 
不 稳定 状态 转变 为 稳 态 则 不 一 定 会 发 生 ， 这 是 由 于 
清 后 现象 的 存在 。 

燃烧 不 稳定 中 出 现 的 分 爹 ， 滞 后 ， 极 限 环 和 触 
发 等 现象 均 可 以 用 动力 系统 理论 来 解释 和 描述 ， 这 
也 是 目前 研究 的 一 大 新 兴 热 点 。 通 过 非 线 性 动力 学 
理论 ， 可 对 非 线 性 系统 的 特征 进行 深入 的 刻画 、 分 
析 和 诊断 。 现 有 研究 大 部 分 以 燃气 轮机 和 火箭 发 动 
机 出 现 的 热 声 不 稳定 问题 为 背景 ， 研 究 对 象 则 通常 
为 各 种 简化 的 燃烧 器 或 实际 的 燃气 轮机 和 发 动机 。 
其 中 , Rijke 管 是 研究 热 声 现象 最 方便 最 典型 的 实验 
系统 ， 也 是 研究 燃烧 室 火焰 区 燃烧 与 声场 的 相互 作 


变 为 大 幅度 自持 振荡 ， 尽 管 系统 是 线性 稳定 的 。 而 
究 证 明了 某 些 类 型 的 噪声 更 能 触发 振荡 。 

现 有 研究 对 Rijke 管 热 声 不 稳定 非 线 性 动力 学 
有 一 定 认 识 ,但 对 系统 分 命 性 质 的 认识 尚 不 够 清晰 ， 
对 于 无 量 纲 加 热 功率 、 加 热 器 位 置 、 阻 尼 系 数 等 系 
统 参数 对 分 命 的 影响 规律 缺乏 了 解 。 本 文 建立 了 
Rijke 管 热 声 不 稳定 的 简化 模型 , 得 到 其 一 维 声场 的 
控制 方程 ， 并 对 其 进行 了 求解 ， 通 过 对 Rijke 管 的 
热 声 不 稳定 性 进行 分 命 分 析 ， 加 深 对 Rijke 管 热 声 
不 稳定 产生 机 理 和 振荡 现象 的 认识 。 
1 热 声 不 稳定 模型 及 求解 
1.1 Rijke 管 模型 

本 文 所 研究 的 水 平 Rijke 管 模型 如 图 1 所 示 。 
在 该 模型 中 ， 管 长 工 是 固定 不 变 的 ， 加 热 器 的 位 置 
好 的 值 可 以 通过 改变 加 热 器 在 管道 内 的 位 置 来 改变 ， 
热源 为 电热 丝 加 热 装置 ， 在 实验 条 件 下 其 加 热 强 度 
可 以 通过 改变 通过 的 电流 来 改变 。 


Heat source 


图 1 水 平 Rijke 管 的 结构 示意 图 ， 来 流 平均 速度 为 x ， 管 长 
为 L， 热 源 与 来 流 开口 端 距离 为 zy。 


Fig.1 Structure schematic of a horizontal Rijke tube. The mean 


flow velocity is u , the length of the duct is L and the distance 


between the heat source and the incoming flow open end is 3. 
1.2 控制 方程 

忽略 平均 流量 与 温度 梯度 的 影响 ， 控 制 管内 声 
场 的 动量 和 能 量 方程 如 下 器: 
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为 了 简化 和 方便 分 析 ， 对 式 (1) 和 (2) 进行 
无 量 纲 化 处 理 ， 得 到 无 量 纲 化 的 系统 控制 方程 组 : 
a a 
“6r Ox 
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A O) - (5) 中 ，x 是 沿 轴 向 方向 的 距离 ，t 是 时 
H p 为 管内 的 平均 压强 ，u' 为 声波 速度 ，p' 为 
声波 压强 ， 2 为 流体 密度 ，7 是 介质 的 比 热 比 ， 声 
速 为 co MM, 为 平均 流 马 赫 数 。 此 外 ， 为 阻尼 系 
Bl, Get) 为 单位 面积 上 由 于 电 加 热 器 产生 的 热 释放 
速率 波动 ，6 为 标准 狄 拉 克 分 布 ，xy 为 热源 位 置 。 
符号 “” 和 ““ ”分别 表示 有 量 纲 量 和 平均 旱 。 
阻尼 可 以 用 下 式 表示 : 
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参数 c M c, 保持 常数 不 变 , w; = jz 为 第 j 阶 
模 态 的 无 量 纲 频率 。 
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热源 模型 采 / 
Ot) = 


aly = x Jace dy x| | 性: uy (t —7) 
SV3cp 2 3 
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其 中 , ,和 4 代表 电热 丝 的 长 度 和 直径 , TN 
电热 丝 的 温度 ，T 代表 周围 空气 的 平均 温度 ，S 为 
管道 截面 积 ，4 代表 空气 的 热传导 率 ，c, 表示 恒定 
体积 内 单位 质量 空气 的 比热容 ， 广 为 气体 的 平均 密 
度 。 由 于 热 惯 性 的 存在 ， 传 热 和 流动 速度 之 间 有 一 
个 时 间 延 迟 z 。 需 要 注意 的 是 ， 此 热源 模型 只 适用 
于 “紧凑 型 ”热源 的 描述 ， 即 热源 厚度 相对 波长 而 
言 可 以 忽略 不 计 。 

将 式 (7) 代入 
可 改写 为 : 


J Heckl 经 验 公式 [19; 


能 量 方程 (4) 中 ， 则 能 量 方程 
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1.3 Galerkin 展开 


求解 Rijke 管 装 置 热 声 问题 的 常规 方法 是 使 用 
Galerkin 投影 和 模 态 扩展 .Galerkin 方法 的 原理 是 任 
何 域 中 的 函数 都 可 以 表达 成 该 域内 展开 函数 的 有 登 加 。 
基础 函数 的 选取 原则 是 它们 须 满足 边界 条 件 ， 且 基 
础 函数 的 选择 不 是 唯一 的 。 在 这 样 的 方法 中 ， 从 满 
足 边 界 条 件 的 空间 基 函 数 来 看 ， 空 间 和 时 间 上 变化 
的 压力 和 速度 信号 被 扩展 ， 这 些 基 函 数 的 选择 不 是 
唯一 的 。 本 文选 择 的 基础 函数 是 任意 的 ， 并 非 系统 
的 特征 函数 ， 而 是 线性 化 后 的 系统 自 共 力 


(self-adjoint) 部 分 的 特征 函数 。 因 此 ， 速 度 场 和 
压力 场 可 以 表示 成 管道 回 有 模 态 的 形式 : 
N 
u' = 》Ui cos(k;x) (9) 
j=l 
N 
p'=7M > P; sin(k jx) (10) 
jal 


A, k= jr 为 第 j 个 管道 模 态 的 无 量 纲 波 长 。 
在 极限 Nw 条 件 下 , 这 些 基 函数 形成 一 个 完整 的 


基础 


ER, 显然 ,本 文选 取 的 展 了 


于 函数 满足 边界 条 件 。 
将 式 (9) Al C10) 替换 到 式 G) 和 8) 中， 
并 投影 (projecting along) 到 基础 函数 系 ， 得 到 : 
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环 的 振幅 和 相位 差别 较 大 ， 但 随 着 Galerkin 模 态 数 
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ASPET, BER STIs A 
据 图 2, 我 们 可 以 看 


对 式 (15) 和 (16) 的 数值 求 , 对 不 同 数目 的 Galerkin 
模 态 下 的 系统 演化 
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[0 相位 差别 逐渐 减 小 。 根 
Hs H O 阶 和 10 阶 声学 模 态 
得 到 数值 解 的 差异 是 几乎 可 以 忽略 不 计 ， 也 再 次 说 
明了 10 阶 声学 模 态 Galerkin 投影 的 收敛 性 。 对 于 
10 阶 声学 模 态 , 增加 一 阶 Galerkin 模 态 所 改变 的 声 
波 速度 极限 环 振幅 小 于 1.4%”。 因 此 , 在 本 文 以 后 
-Gan 的 所 有 计算 中 ,我 们 采用 10 阶 Galerkin 模 态 的 模型 
来 进行 线性 和 非 线 性 稳定 性 分 析 以 确保 收敛 性 。 
(12) 
其 中 ， 无 量 纲 加 热 功率 为 : = 
1-0 F —e— 9 Galerkin modes j 
_ Ay -Db, T -Py fopa ae || —o~ 10 Gelaerkin mode 
yMaScy DA 
上 述 完 全 耦合 的 Galerkin 模型 可 以 通过 线性 化 00 
延迟 进一步 简化 ， 该 方法 只 适用 于 满足 r<7 的 
Galerkin KEF, FT, = 2/7 为 第 j 阶 Galerkin 模 
态 的 时 间 周 期 。 在 满足 这 种 假设 的 条 件 下 ， 延 迟 项 
© 可 以 写成 : Time ran (0) 7 
— U,(t—t) SU, +tk,P. (14) 
N 通过 线性 化 延迟 ， 得 到 最 终 Galerkin 展开 的 常 图 2 不 同 声学 Galerkin 模 态 数目 下 系统 随时 间 演化 的 比较 ， 
v= 微分 方程 组 ; 系统 参数 为 : K=1,c1=0.1, 00.06, x=0.3, 1=0.2. 
= Fig.2 Comparison of time evolution of the system with different 
a t oA = _ acoustic Galerkin mode numbers. The parameter values of the 
B+ 26 ,@;P; 一 kU, system are K=1, cl=0.1, c2=0.06, x=0.3 and t = 0.2. 
= ‘| | LS east +k,P)| - sina yxp) 2 线性 稳定 性 分 析 
= 为 研究 无 穷 小 扰动 对 系统 稳定 性 的 影响 ， 首 先 
Gabon) (16) ”对 系统 进行 线性 稳定 性 分 析 。 如 果 系统 逐渐 远离 
定 状态 ， 则 系统 最 终 是 不 稳定 的 ， 如 果 系 统 逐 渐 接 
其 中 ，k; = jz 为 第 j 阶 管内 声 模 态 的 无 量 纲 波 数 。 近 稳 定 状 态 
1.4 Galerkin 模 态 收敛 
在 极限 N 一 条件 下 ， 式 〈9) 和 C10) 这些 
基 函 数 形成 一 个 完整 的 基础 系 ， 但 是 考虑 到 计算 的 
可 行 性 ， 我 们 只 能 使 用 有 限 数 量 的 模 态 数目 。 精 确 
捕捉 到 系统 的 线性 丰 


局 部 分 析 系 统 在 稳定 状态 附近 时 
一 方 


ea 
I 


1 非 线性 行为 所 需 的 Galerkin 模 
态 的 数目 被 称 为 模 态 收敛 。 利用 Matlab 编程 实现 了 


ERE PES BT EE 


的 稳定 性 变化 。 另 


行 比较 。 对 于 在 线性 不 稳定 区 


究 有 限 振幅 扰动 对 
响 效果 ， 并 用 于 表征 所 得 到 的 系统 渐 近 状 
临界 点 称 为 分 岔 点 。 一 


域 ， 具 有 不 同 数量 声学 模 态 的 该 系统 的 时 间 演 化 的 
比较 见 图 2。 可 以 得 到 1 阶 Galerkin 声学 模 态 下 极 


被 确定 ， 该 分 岔 点 本 身 继续 与 系统 另 一 个 相关 参数 
变化 。 分 岔 点 分 支 提 


态 。 由 于 系统 平衡 解 失 稳 ， 系 统 的 行为 发 生 改 变 的 


限 环 的 振幅 和 相位 与 10 阶 Galerkin 声学 模 态 下 极限 


日 一 个 分 岔 点 根据 茶 个 参数 
eT eet 
区 域 。 这 个 

一 个 


可 以 分 离 相关 参数 空间 中 线性 稳定 和 线性 不 稳定 的 
i 


稳定 的 边界 ， 该 边界 
当 的 二 维 投影 上 


稳定 边界 是 在 所 有 自由 参数 变化 空间 的 
超 曲面 ， 但 比较 方便 的 是 我 们 可 以 表示 几 个 适 
Hk- 
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在 Rijke 管 中 ， 系 统 参数 有 加 热 器 的 无 量 纲 加 
热 功 率 、 加 热 器 的 位 置 、 阻 尼 系 数 以 及 时 间 延 迟 ， 
但 是 在 实际 实验 操作 条 件 下 ， 可 以 精确 改变 的 参数 
只 有 无 量 纲 加 热 功率 和 加 热 器 的 具体 位 置 ， 时 间 延 
述 和 阻尼 系数 都 是 根据 其 他 条 件 来 改变 的 。 我 们 以 
加 热 器 位 置 和 时 间 延 迟 项 的 线性 稳定 性 边界 为 例 ， 
对 加 热 器 位 置 和 时 间 延 迟 变化 时 的 线性 稳定 性 区 域 
和 线性 不 稳定 区 域 进行 分 析 。 一 个 典型 的 加 热 器 的 
位 置 和 时 间 延 迟 的 线性 稳定 边界 变化 如 图 3 所 示 ， 
显示 了 对 于 确定 的 阻尼 和 加 热 器 功率 等 固定 系统 参 
数值 ， 对 于 加 热 器 位 置 (xy) 的 一 个 选择 范围 ， 系 


TT 
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型 图 4 所 示 。 如 果 极 限 环 是 不 稳定 ， 如 图 4 Ca) 所 
示 ， 不 稳定 极限 环 分 支 形成 一 个 双 稳 态 区 域 ， 在 该 
区 域内 稳 态 解 对 于 小 扰动 是 稳定 的 而 对 于 大 扰动 是 
不 稳定 的 。 然 而 ， 不 稳定 极限 环 的 这 一 分 支 可 能 会 
BREE? AME Ee. MPD RNB 
值 ， 稳 定 解 对 于 任何 幅度 的 扰动 都 是 稳定 的 ， 因 此 
系统 是 全 局 稳定 的 。 因 此 在 亚 临 界 分 兮 情况 下 ， 线 
性 (局部) 和 非 线 性 (全 局 ) 的 稳定 性 边界 是 不 同 
的 。 这 种 分 命 行为 被 称 为 亚 临 界 Hopf 分 贫 。 但 如 
果 极 限 环 是 稳定 的 如 图 4 (b) 所 示 ， 系 统 平稳 地 从 
稳定 的 稳 态 解 发 展 为 一 种 具有 递增 极限 环 振幅 不 稳 


统 是 否 为 线性 不 稳定 则 取决 于 时 间 延 迟 Cr )， 反 之 
亦 然 。 对 于 较 小 或 合理 大 的 rz ， 比 如 r<0.15 和 
T>0.85, ， 加 热 器 在 任何 位 置 时 ，Rijke 管 系统 都 是 
线性 稳定 的 ， 即 对 于 小 扰动 ， 系 统 都 最 终 趋 于 稳定 
状态 。 只 有 在 0.15<r< 0.85 范围 内 ， 系 统 平衡 解 的 
稳定 性 取决 于 加 热 器 的 位 置 。 


— 


f Linearly stable region 


0.0 0.1 0.2 0.3 0.4 0.5 


图 3 加 热 器 位 置 六 和 时 间 延 迟 r 的 线性 稳定 性 边界 , 系统 其 
也 参数 为 : K=0.8，c1=0.1，c2=0.06。 


Fig.3 Linear stability boundary between xf and t. The other 


parameter values of the system are K=0.8, ci=0.1 and c2=0.06. 
3 分 义 分 析 

DEAT EH FED BINA ATE EE APE LEE 
分 岔 现象 是 指 对 于 某 些 完全 确定 的 非 线性 系统 ， 当 
RANK BR u 连续 变化 到 茶 个 临界 值 凡 时 , 系 
统 的 全 局 性 性 态 〈 定 性 性 质 、 拓 扑 性 质 等 ) 会 发 生 
突然 变化 。 公 称 为 参数 u 的 分 贫 值 或 分 支 值 。 系 统 
稳定 状态 失去 稳定 性 出 现 孤 立 的 极限 环 周 期 解 的 分 
BAA Hopf 分 倪 。 极 限 环 的 新 兴 分 支 的 稳定 性 决 
ET Hopf 分 支 的 性 质 或 类 型 。 这 两 种 Hopf IAK 


定 的 稳 态 解 。 这 种 分 岔 类 型 称 为 超 临 界 Hopf 分 贫 。 


(b) 2 


Measure 


0.5 
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图 4 (a) 亚 临界 和 Cb) 超 临界 Hopf 分 岔 附近 参数 变化 时 
对 应 的 分 盆 行 为 n31。 


Fig.4 Bifurcation behavior of a measure for the variation of a 


parameter near (a) Subcritical and (b) Super-critical Hopf 
bifurcation. 

可 通过 定性 分 析 方 法 判定 亚 临 界 和 超 临 界 
Hopf DÉ (1+ X)* 形式 的 源 项 非 线性 分 岔 的 性 质 
可 以 通过 把 非 线性 项 X 展开 成 级 数 , 并 放弃 高 阶 项 
来 判断 。 上 述 (1+X) 表达 式 的 二 项 扩展 的 结果 为 


下 面 的 式 子 : 
jeep OD yoy Ole?) ys 


(17) 

在 上 述 表达 式 中 ， 当 0<|a|<1 和 |a|>2 时 ,一 阶 项 
aX 和 三 阶 项 a(g 一 1)(ga 一 2)X /3! 的 符号 是 相同 
的 ， 而 当 1<|a|<2 时 ， 一 阶 项 aX 和 三 阶 项 
a(o-1)( Q-2)X3/3! 这 两 项 的 符号 是 不 同 的 ,一 阶 
和 三 阶 项 的 符号 标志 决定 了 分 命 点 的 性 质 。 每 当 这 
两 项 具有 相同 的 符号 时 , 分 岔 是 亚 临界 分 岔 类 型 的 ， 
而 当 这 两 项 具有 不 同 的 符号 时 ， 分 侈 为 超 临界 分 贫 
类 型 的 。 对 于 Rijke 管 模型 中 的 热 释 放 速 率 波动 ， 
Q =1/2， 这 意味 着 该 模型 将 出 现 亚 临界 类 型 Hopf 
NB o 

在 Rijke 管 热 声 系统 中 ， 可 变 的 系统 参数 有 有 
加 热 器 的 无 量 纲 加 热 功 率 K、 加 热 器 的 位 置 x BE 
尼 系 数 c 以 及 时 间 延 迟 z 。 分 别 对 系统 内 无 量 纲 加 
热 功率 及 、 加 热 器 的 位 置 x 阻尼 系数 c 以 及 时 间 
延迟 zt 的 变化 进行 分 岔 分 析 。 
3.1 加 热 器 功率 的 影响 

E Rijke 管 热 声 系 统 中 的 加 热 器 为 电热 丝 加 热 ， 
通过 改变 K 值 来 探究 加 热 器 功率 对 系统 产生 的 影 
响 。 无 量 纲 加 热 器 功率 可 以 通过 增加 提供 给 加 热 器 
的 功率 来 达到 ， 其 增 大 代表 着 系统 驱动 力 的 增加 ， 
增加 的 驱动 力 使 得 系统 变 得 更 加 不 稳定 。 因 此 ， 对 
于 较 小 的 天 值 ， 系 统 的 平衡 是 稳定 的 ， 在 所 有 初始 
扰动 下 ， 都 会 衰减 渐 近 接近 于 零 。 增 加 天 值 将 降低 
流动 的 稳定 性 裕 度 ， 在 分 命 点 附近 ， 系 统 演变 为 线 
性 不 稳定 ， 从 而 导致 管内 的 振荡 流 模式 。 图 5 Nit 
算得 到 的 无 量 纲 加 热 器 功率 (KK) 的 变化 对 系统 演 
化 的 影响 。 空 心 圆 表明 不 稳定 解 ， 实 心 圆 代表 稳定 
解 。 由 图 可 知 靠近 Hopf 分 命 点 的 小 振幅 极限 环 是 
不 稳定 的 ,这 些 不 稳定 的 极限 环 与 稳定 平衡 解 共存 。 
这 些 极限 环 的 不 稳定 分 支 进 一 步 经 历 一 个 折合 点 分 
命 获 得 稳定 。 据 此 可 判定 分 兮 属于 亚 临 界 类 型 ， 证 
实 了 之 前 的 结论 。 在 其 他 系统 参数 一 定 的 条 件 下 ， 
系统 随时 间 的 演化 行为 取决 于 K 值 的 变化 。 在 
<0.515 时 ,系统 对 于 任何 幅度 的 扰动 ， 都 最 终 趋 
于 零 ， 此 区 域 为 全 局 稳定 区 域 。 随 着 天 值 的 不 断 增 
大 ， 当 0.515< 玉 < 0.615 时 ， 系 统 有 三 种 可 能 的 状 
态 ， 稳 定 状 态 ， 不 稳定 极限 环 以 及 稳定 极限 环 。 根 
据 初 始 条 件 的 不 同 ， 系 统 可 能 最 终 进入 稳定 状态 或 
者 周期 振荡 ， 该 区 域 被 称 为 双 稳 态 区 域 。 当 
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玉 >0.615 时 ， 系 统 进 入 全 局 不 稳定 区 域 , 对 于 任意 
小 幅度 的 初始 扰动 ,系统 最 终 都 会 达到 极限 环 状态 。 
通过 无 量 纲 加 热 功率 的 分 岔 图 也 可 以 证 明 无 量 纲 加 
热 功率 的 增 大 是 系统 产生 热 声 不 稳定 的 主要 原因 。 


O unstable solutions 
@ stable solutions 
0.8 È 
e 
e 
0.6 ee 
= e 
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Fold bifurcation O Hopf bifurcation 
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图 5 无 量 纲 加 热 器 功率 天 变化 时 的 分 岔 图 ， 系 统 其 他 参数 
值 : cl=0.1，cz=0.06，x=0.3，r= 0.2。 


Fig.5 Bifurcation plot for variation of non-dimensional heater 
power K. The other parameter values of the system are ci=0.1, 
c2=0.06, x¢=0.3 and t = 0.2. 
3.2 阻尼 系数 的 影响 

为 了 研究 系统 中 阻尼 变化 对 系统 响应 的 影响 ， 
阻尼 模型 中 的 一 个 阻尼 参数 〈c; ) 是 可 以 变化 的 。 
在 实际 的 实验 条 件 下 ， 系 统 阻 尼 的 大 小 是 通过 改变 
管道 的 末端 条 件 来 改变 。 根 据 非 线性 理论 ， 增 加 阻 
尼 可 以 增强 系统 的 稳定 性 。 图 6 为 系统 行为 演化 随 
阻尼 系数 (ci ) 的 变化 的 分 贫 图 谱 。 如 预期 的 一 样 ， 
增加 阻尼 确 具 有 增强 系统 稳定 的 作用 ， 因 为 平衡 解 
对 于 较 大 的 阻尼 系数 ， 取 任意 时 间 延 迟 z 都 是 稳定 
的 ， 降 低 阻 尼 则 使 系统 失去 稳定 性 。 对 于 图 示 的 时 
间 延 迟 值 、 无 量 纲 加 热 功 率 以 及 加 热 咒 位 置 参 数 ， 
存在 一 个 临界 值 = 0.19 ， 当 低 于 该 值 时 ， 系 统 对 
于 任何 幅度 的 初始 扰动 最 终 都 会 发 展 成 稳定 振幅 的 
极限 环 ， 该 临界 值 即 为 亚 临 界 Hopf PA. 
c <0.19 时 ， 对 于 任何 幅度 的 初始 扰动 ， 系 统 最 终 
都 会 发 展 成 具有 稳定 振幅 的 极限 环 ; 当 
0.19<c <0.25 时 ， 该 区 域 为 双 稳 态 区 域 ， 在 该 区 
域内 ， 对 于 较 大 幅度 的 初始 扰动 ， 系 统 发 展 成 为 稳 
定 极限 环 ， 而 对 于 小 振幅 的 扰动 ， 系 统 最 终 会 趋 于 
平衡 ; 随 着 阻尼 系数 ci 的 进一步 增 大 , 当 c > 0.25 时 ， 
对 于 任意 幅度 的 初始 扰动 ， 系 统 都 最 终 趋 于 稳定 状 
S> c =0.25 AHER. 
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6 阻尼 系数 (ct) 变 化 时 的 分 岔 图 , 系统 的 其 
天 =1，c2=0.060，Xf=0.3，r= 0.2。 


它 参 数值 为 : 


Fig.6 Bifurcation plot for variation of damping coefficient c1. 
The other parameter values of the system are K=1, c2=0.06, xf 
=0.3and t= 0.2. 
3.3 加 热 器 位 置 的 影响 

加 热 器 的 位 置 Op) 对 系统 动力 学 也 有 显著 的 
影响 ， 加 热 器 位 置 的 改变 是 通过 把 加 热 器 放置 在 治 
管道 长 度 不 同 的 位 置 处 来 实现 。 系 统 稳定 性 的 变化 


以 一 种 特别 的 方式 取决 于 管道 内 加 热 器 的 位 置 变化 。 


加 热 器 位 置 (xy) 变化 时 的 分 岔 图 谱 如 图 7 所 示 。 
当 管道 内 加 热 器 的 位 置 从 上 游 的 开端 开始 变化 时 ， 
该 系统 最 初 是 线性 稳定 的 。 在 加 热 器 的 位 置 zz 的 临 
界 值 ， 系 统 演变 为 线性 不 稳定 。 随 着 加 热 器 位 置 的 
步 变化 ， 系 统 仍然 是 线性 不 稳定 的 ， 直 到 加 热 
器 位 置 位 于 xp， 在 该 点 形成 男 一 个 Hopf TA R, 
系统 也 重新 变 成 稳定 状态 。 加 热 器 位 置 Cp 变化 
的 分 岔 图 表明 ,在 管 的 两 端 ，Hopf DE AAE 
F Hopf 分 岔 。 在 两 个 Hopf 分 岔 点 之 间 ， 存 在 一 个 
全 局 不 稳定 区 域 。 在 全 局 不 稳定 的 区 域内 ， 对 于 任 
何 的 初始 条 件 ， 系 统 将 渐 近 达到 相应 稳定 振幅 的 极 
限 环 。 根 据 加 热 器 位 置 变化 时 的 分 岔 图 ， 我 们 可 以 
得 到 两 个 亚 临界 Hopf 分 兮 点 的 具体 值 。 对 于 我 人 
选 定 的 其 他 系统 参数 ， 当 xj <0.11 时， 对 于 任何 幅 
度 的 初始 扰动 ， 系 统 最 终 都 会 趋 于 稳定 ; 4 
0.11<xr <0.14 时 ， 系 统 处 于 双 稳 态 区 域 ， 系 统 根 
据 初 始 条 件 的 不 同 产生 不 同 的 演化 行为 ， 则 第 一 个 
亚 临 界 Hopf D a AME xn=0.1; 当 
0 . I< 和 rj < OR, 系统 处 于 全 局 不 稳定 区 域 , 对 
于 任何 的 初始 条 件 ， 系 统 将 渐 近 达到 相应 稳定 振幅 
的 极限 环 ; 当 0.365 < x, <0.39 时 ,系统 又 进入 双 稳 
态 区 域 ， 则 第 二 个 亚 临 界 Hopf 分 岔 点 位 置 
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xj) =0.36;， 随 着 加 热 器 位 置 进一步 变化 ， 当 
x, > 0.39 时 ， 系 统 重新 进入 全 局 稳定 性 区 域 ,对 于 
任何 幅度 的 初始 扰动 ， 系 统 最 终 都 会 趋 于 稳定 。 
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7 加 热 器 位 置 (xy) 变 化 时 的 分 命 图 , 系统 的 其 它 参 数 为 : 
c1=0.1, c2=0.06, K=0.8, T=0.2。 


Fig.7 Bifurcation plot for variation of location of heater (x). 
The other parameter values of the system are cl = 0.1, c2 = 0.06, 
K=0.8 and t=0.2. 
3.4 时 间 延 迟 的 影响 

在 Rijke 管 中 ， 由 于 加 热 器 与 流动 介质 之 间 热 
惯性 (thermal inertia) 的 存在 ， 传 热 和 流动 速率 的 
变化 之 间 存 在 一 个 时 间 延 迟 z ， 其 大 小 与 加 热 器 中 
电热 丝 的 直径 以 及 来 流速 度 有 关 。 在 实际 实验 条 件 
Po WEIER Cr) 很 难 精确 的 进行 改变 。 时 间 
HEIR Cr) 变化 时 对 应 的 分 贫 图 谱 见 图 8。 亚 临界 
Hopf 分 岔 点 对 应 的 时 间 延 迟 值 为 rz=0.155。 当 
z<0.13 时 ， 系 统 为 处 于 线性 稳定 状态 ， 对 于 任意 
晶 度 的 初始 扰动 ， 系 统 最 终 都 趋 于 稳定 状态 ; 当 
0.13<z<0.155 时 ， 系 统 对 于 小 振幅 的 扰动 最 终 趋 
于 平衡 ， 而 对 于 大 振幅 的 扰动 ， 系 统 发 展 成 为 极限 
环 ， 即 双 稳 态 区 域 ， 当 +>0.155 时 ， 系 统 处 于 线性 
不 稳定 区 域 ， 对 于 任意 小 的 初始 扰动 ， 系 统 最 终 都 
会 发 展 成 为 极限 环 。 值 得 注意 的 是 ， 发 现在 
5>0.155 的 线性 不 稳定 区 域内 , 随 着 时 间 延 迟 Cr ) 
MEHAR, PRIE |U, | 的 值 先 逐 渐 增 大 后 又 逐渐 减 
小 。 这 一 现象 可 通过 瑞 利 准则 来 解释 07。 时 间 延 迟 
的 存在 是 因为 热源 的 放 热 和 流动 之 间 存 在 热 惯性 ， 
当 加 热 器 功率 增 大 时 ， 时 间 延 迟 会 随 之 增 大 。 时 间 
延迟 的 增 大 使 得 声场 振荡 和 放 热 脉动 的 相位 差 发 生 
变化 ， 当 二 者 相位 差 逐渐 增 大 而 远离 0” 时 ， 振 荡 
局 度 会 逐渐 减 小 ， 从 而 形成 了 这 一 现象 。 据 此 也 可 
以 推测 当 z=0.4 时 ， 声 场 振荡 和 放 热 脉动 的 相位 


rl 


= 


差 刚好 为 0”， 故 此 时 振荡 幅 值 最 大 。 


4 结论 


本 文 对 Rijke 管 热 声 系 统 进行 建 模 ， 利 用 
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8 时 间 延 迟 (z) 变 化 的 分 贫 图 , 系统 的 其 它 参数 为 : c1=0.1， 
c2=0.06, K=0.8, x=0.3. 


Fig.8 Bifurcation plot for variation of time lag (z). The other 
parameter values of the system are ci=0.1, c2=0.06, K=0.8, xf 


=0.3. 
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